DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


Development  of  Turbulent  Biological  Closure  Parameterizations 

Louis  Goodman 

School  for  Marine  Science  and  Technology  (SMAST) 

University  of  Massachusetts  Dartmouth 
706  South  Rodney  French  Blvd 
New  Bedford,  MA  02744 

phone:  (508)  910-6375  fax:  (508)  910-6376  email:  lgoodman@umassd.edu 

Award  Number:  N000 1408 10631 
http://www.smast.umassd.edu/Turbulence/ 


LONG-TERM  GOAL 

The  long-term  goals  of  this  project  are: 

(1)  to  develop  a  theoretical  framework  to  quantify  turbulence  induced  NPZ  interactions. 

(2)  to  apply  the  theory  to  develop  parameterizations  to  be  used  in  realistic  environmental 
physical  biological  coupling  numerical  models. 

OBJECTIVES 

Connect  the  Goodman  and  Robinson  (2008)  statistically  based  pdf  theory  to  Advection  Diffusion 
Reaction  (ADR)  modeling  of  NPZ  interaction. 

APPROACH 

A  nonlinear  model  for  biological  and  physical  dynamical  interactions  in  a  laminar  flow  field  being 
upwelled  into  the  mixed  layer  Robinson  (1999)  (Fig  1)  has  been  extended  to  turbulent  flow  (Goodman 
and  Robinson,  2007).  The  approach  of  the  Goodman  and  Robinson  theory  has  been  to  develop  a 
probability  density  function  (pdf)  for  the  turbulent  displacement  field  and  use  that  to  calculate  the 
turbulence  induced  biological  interaction  (TIBI)  terms,  i.e.  <  PNJ  >,<  PZ ,  >,<  ZiNj  >,  where 

N  i,Pi,Zi  are  the  ith  component  of  a  field  of  different  nutrients,  phytoplankton,  and  zooplankton 
embedded  in  the  turbulent  field.  Contrast  the  TIBI  terms  with  the  biological  turbulent  flux  terms 
<u'Ni>,<u'Pi>,<u'Zi>.  The  formalism  for  modeling  the  latter  type  of  terms  are  well  developed, 
typically  involving  some  type  of  eddy  diffusivity  or  higher  order  closure  such  as  Mellor  and  Yamada, 
(1982).  However,  at  present,  no  biodynamical  basis  for  closure  of  the  TIBI  terms  has  been  developed. 

The  current  approach  to  handling  the  TIBI  terms  is  either:  (1)  to  ignore  them  by  setting 
<  (••),(••),  >=<  (••),  ><  (••),  >  in  an  advection  diffusion  reaction  (ADR)  equation  approach  (Donaghay, 

and  Osbom,  1997),  or  (2)  to  perform  a  numerical  simulation  for  the  turbulent  displacement  field  and 
explicitly  calculate  the  TIBI  terms.  The  former,  as  we  have  shown  (Goodman  and  Robinson,  2007), 
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can  result  in  a  large  overestimate  of  the  effect  of  turbulence  on  these  interactions,  while  the  latter  is 
very  limited  in  its  domain  size,  suffers  from  difficulty  in  imposing  realistic  boundary  condition  at  the 
laminar  turbulent  interface,  and  can  only  reveal  significant  physics  of  the  TIBI  terms  with  a  large 
number  of  repetitive  runs,  which  prohibit  the  size  of  the  computational  domain  to  very  limited 
environmental  scenarios.  In  addition,  turbulence  numerical  models  such  as  LES  and  DNS  are  difficult 
to  embed  in  larger  regional  scale  biophysical  coupling  models.  What  is  needed  to  be  useful  in  the 
larger  scale  physical  biological  coupling  models  is  development  of  realistic  parameterizations  of  the 
TIBI  terms,  analogously  to  the  development  of  turbulent  flux  parameterizations  used  in  regional  and 
large  scale  ocean  circulation  models.  In  addition  development  of  such  parameterizations  will  lead  to 
new  physical/  biological  insights  into  the  role  of  the  TIBI  terms. 


Turbulence 


Transition  Zone 


Optically  & 
Biologically  Active 


Laminar  Flow 


x 


> 


Figure  1.  Upwelling  flow  field  into  an  optically  and  biologically  active  mixed  layer  used  in  the 
Goodman  and  Robinson  (2007)  biodynamical  model.  Note  as  a  increases  or  as  the  uptake  time 

increases  relative  to  the  advection  time 


WORK  COMPLETED 

We  have  made  a  major  breakthrough  and  have  obtained  an  exact  solution  to  the  ADR  equation  for 
realistic  boundary  conditions  using  the  Goodman  and  Robinson  pdf  theory  approach.  This  work  is 
now  published  in  the  recent  Issue  of  Dynamics  of  Atmospheres  and  Oceans  dedictaed  to  Allan 
Robinson  who  passed  away  in  2010.  (Goodman  ,L.,  2011) 

To  illustrate  an  application  of  the  theory  and  the  role  of  TIBI  tenn  in  the  ADR  equation  a  simple 
example,  that  of  upwelling  of  seed  nutrients  and  phytoplankton  into  a  turbulent  optically  active  mixed 
layer  is  presented. 

Consider  the  simple  example  used  in  previous  manuscripts  by  Robinson  (1999)  and  Goodman  and 
Robinson(  2008),  namely,  that  of  a  linear  strain  upwelling  of  seed  nutrients  and  phytoplankton  into  an 
optically  active  turbulent  mixed  layer.  The  equations  for  this  model  are 
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where  P  ,  N  are  the  phytoplankton,  nutrient  linear  densities  (units  of  m'3),  respectively;  p  is  the 
nutrient  uptake  rate;  w  =  -az  is  the  vertical  velocity  with  a  constant  strain  rate  ,-a;  z  is  taken  as 
positive  downward.  The  turbulent  mixed  layer  is  located  between  z  =  0  and  z  =  D.  The  mean 
flow  field  is  two  dimensional  and  incompressible.  It  is  assumed  that  one  dimensional  turbulent 
mixing  dominates,  i.e.  turbulent  scales  in  the  horizontal  are  much  larger  than  that  in  the  vertical. 
Adding  (14a)  and  (14b)  result  in  the  total  biomass  density  being  conserved  in  a  Largangian 
coordinate  system,  i.e. 


M0=N  +  P, 


(2) 


where  M0  =  N0+P0  is  the  biomass  density  at  t  =  t0 .  Nonnalizing  the  variables  in  equations  (la) 
and  (lb)  using  the  mixed  layer  depth  D,  total  biomass  density,  M0 ,  and  the  nutrient  uptake  time 
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and  results  in  the  normalized  set  of  equations 

dP  dP  dP 

—  =  ( — )  -az( — =  PN 
dt  dt  ~  dz  1 

dN  ,57V  dN , 

- =  ( - )7  -  az( - ).  =  -  PN, 

dt  dt  ~  dz  1 


D 


which  have  the  solutions 


P  =  P(Z,t,t„) 
N  =  l-P, 


P0+N0exp[-(t-t0)]’ 


(3a) 

(3b) 


(4a) 

(4b) 


with  P0,N0  the  normalized  phytoplankton  and  nutrient  seed  densities,  respectively,  at  t  =  t0. 
Consider  the  PDF  equation  for  a  linear  strain  rate  mean  flow  and  a  constant  vertical  eddy 
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diffusivity,  k.  .  The  advection  diffusion  (AD)  equation  which  we  assume  the  turbulence 
displacement  probability  density  function  p  satisfies  is  given  by 


/-S  /V  /-\  ZV  /-S  9  ^ 

dp  dp  d  p  _ 

-J—-az-!—-K — V  =  0. 

dt  dz  dz~ 

(5) 

For  the  mixed  layer  boundary  conditions  we  assume  that  the  total  flux  of  material  (the  sum  of  the 
advective  and  turbulent  components)  vanish  at  the  base  of  the  mixed  layer  and  both  the  advective 

and  turbulent  flux  vanish  at  the  top  of  the  mixed  layer,  whence 

/-N  /V  A  A 

dp  dp 

azp  +  K — =  k —  =  0  atz  =  0, 

dz  dz 

(6a) 

dp 

azp  +  K —  =  0  atz  =  l. 

dz 

(6b) 

Equations  (6a)  and  (6b)  insure  that  a  normalized  solution  of  p  can  be  interpreted  as  a  PDF. 
Equation  (5)  with  boundary  conditions  (6a)  and  (6b)  also  describes  the  evolution  of  a  scalar  non 
interacting  density  field,  p  .  To  see  the  latter,  consider  a  solution  to  (5)  of  the  form 


P  =  J  dfK(f)exp(-ft)  +  p0, 

0 


where  p0  is  the  density  at  z  =1,  t  =  to  and  which,  upon  substitution  into  (5),  yields 


dK  82K 

vK  +  az - yk — — 

dz  dz2 


with  the  boundary  conditions  (6)  becoming 


dK 

dz 


=  0  at  z  =  0, 


dK 


Let 


aK  +  k —  =  0  at  z  =  1 . 

dz 

Pe  2 


K  =  Gexp( - z  ) 


where  the  turbulent  Peclet  number  is  given  by 

Pe=°. 

K 


(8a) 

(8b) 
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Substitution  of  K  into  equation  (7)  yields 
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where  Y  —  CCy  ,  With  the  above  substitution,  G  satisfies  the  boundary  conditions 


—  zz  0  at  z  =  0, 

dz 


Pe  ^  dG 
2  dz 


=  0 


at  z  =  1 . 


(10a) 

(10b) 


It  is  straightforward  to  show  that  Eq.  (9)  with  boundary  conditions  (10)  result  in  eigenfunction 
solutions  G  =  Gm  ,  with  associated  eigenvalues,  y  =  ym  ,m  =  1,2 Note  that  G;;;is 
orthononnal, 

l 

[dzG  G  =8  . 

I  m  n  mn 
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It  is  also  straightforward  to  show  that 


g, 
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using  these  eigenfunctions  and  eigenvalues  yields  the  solution 


P  =  exP  (1 1) 

where 


5 


jtfz'G„,(z')exp(^z'2). 


Equation  (11)  describes  the  temporal  and  spatial  evolution  of  the  density  field,  p  ,  which  is  given  by 
p  =  p0  at  t  =  to  and  z  =1 .  Note  that  when  t  =  oo,  p  =  pQ ,  as  expected. 

To  obtain  the  PDF,  we  rewrite  (1 1)  as 

dt'H(to,t%~)=  }  df{p0H(t0,t')}{—(-^)},  (12) 

L  Po  dt 


with  the  Heaviside  function  H  defined  by 


1 

0 


t'>t0 

t'<t 0 


Eq.  (12)  can  then  be  interpreted  as  the  prescription  for  obtaining  the  average  density  P  ~  P  from  the 

initial  density  p  =  p0H(t0,t ')  using  the  PDF,  F  =  -(-—)  . 

Po  dt' 


Changing  the  independent  random  variable  C  to  z  according  to  z=exp (-a(t-t')  ,  results  in 


P  =  J  dzp(z)  =  J  dzp(z)F(z), 

0  _  0 


dt 1 


F  =  ym=y  A  G  cxGxz 

/  j  m-\  J  m  m  m  1  1 


(r,„-i)  _ 


50 


with 


=°° 

^  //7=1  m  mi  1 


Note  that  F,  as  a  PDF,  has  the  proper  normalization 


|  dt  'F(t ')  =  J  dzF(z)  =  jdS^-  =  g(l)  -  g(0)  =  1 

n  0Z 


(13a) 

(13b) 
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We  can  interpret  z  as  the  Lagrangian  position  of  a  fluid  particle  being  upwelled  by  the  linear  strain 
vertical  velocity,  w  =  -dz,  into  the  turbulent  mixed  layer  at  time  t,  whose  initial  position  at  z  =  1  was  f 
.  Ensemble  averages  of  P  and  PN  are  then  given  by 

t  1 

<  P  >=  f  dt  P(t  ')F(t  ’)  =  f  dzP(z)F(z),  (14a) 

—oo  0 

t  1 

<  PN  >=  |  dt  P(t  ')N(t  ')F(t  0=1  dzP(z)N (z)F(z).  (14b) 

—00  0 


Using  equation  (12  b)  with  (14a)  and  (14  b),  it  is  straightforward  to  show  that  the  mean  phytoplankton 
density  P  satisfies  the  ADR  equation 


(/)  (//)  (III)  (IV) 

8P  8P  82P  „ 

- az - k — —  =  R, 

dt  dz  dz~ 

(15)wlth 

boundary  conditions 

dP 

=  0  at  z  =  0, 

dz 

(16a) 

-  dP 

PeP  H - =  PeP  at  z  =  0. 

8z  0 

The 

(16b) 

‘reaction”  term,  R,  is  given  by 

R  =<  PN  >=  Rq+Rt  ; 


with  R0  =  PN  and  R,  =<P'N'>  ,  the  latter  the  TIBI  term.  Thus,  we  have  formally  solved  the  ADR 
equation  (15)  without  recourse  to  dropping  the  TIBI  tenn. 

To  examine  the  role  of  the  TIBI  tenn,  RT  ,  we  will  compare  the  results  out  lined  above  to  that  of 
obtaining  a  solution  to  the  ADR  equation  (15)  with  no  TIBI  tenn,  i.e.  R  =  R0  =  PN  .  To  obtain 
numerical  results  we  use  as  input  parameters  typical  oceanic  values  of 

«  □  10'5 -10'6  sec'1, 

PU  1CT2  -1CT4  — , 
sec 

DU  10-50  m 
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(Large,  1994;  Goodman  and  Robinson,  2008).  This  yields  a  range  of  turbulent  Peclet  numbers 
between  .01  <  Pe  <  10 .  Note  that  Pe  is  independent  of  the  nutrient  uptake  time,  t  .  .  If  we  use  one  day 
as  a  characteristic  nutrient  uptake  time  and  the  larger  value  of  the  linear  strain  rate  of  aU  1 0  5 ,  this 
results  in  the  nonnalized  strain  rate  a  □  .1  (Robinson,  1999). 


In  figure  2  vertical  profiles  of  P  are  shown  for  the  steady  state  (t—>  go, —  =  0  )  for  0  <  Pe  <  go  and 

ct 

a  =  .1 .  The  solid  color  coded  lines  are  Case  A,  the  solution  for  <P>  using  equations  (13b)  and  (14a), 
TIBI  term  included.  The  dashed  color  coded  lines  are  solutions  to  the  ADR  equation  (15)  with  the 
TIBI  term  neglected, RT  =  0  ,  Case  B,  using  the  same  range  of  Pe  and  a -  A.  The  red  solid  line  of 

figure  2  corresponds  to  no  turbulence,  pure  advection,  Pe  =  go  .  We  see  a  very  significant  difference  in 
the  limiting  value  of  Pe  =  0  for  the  two  cases.  Figure  2  shows  that  neglecting  the  TIBI  term  over  a 
wide  range  of  Pe,  up  to  Pe-  10  ,  where  advection  dominates  turbulence,  results  in  a  large  over 
estimate  of  P  .  Compare  the  dashed  with  the  solid  lines.  Thus,  turbulence  in  this  simple  model  tends 
to  limit  phytoplankton  growth  over  a  very  wide  range  of  Peclet  numbers. 

To  understand  this  result  we  present  in  figures  (3a), (3b),  (3c),  and  (3d)  ,  for  the  steady  state  cases 
presented  in  figure  1,  the  contribution  of  the  three  steady  state  terms  of  equation  (15)  -  advection,  term 
II,  blue  lines;  turbulent  diffusion,  term  III,  green  lines;  and  growth  rate  (reaction),  R,  term  IV,  the 
black  lines.  As  in  figure  2,  the  solid  lines  refer  to  Case  A,  the  dashed  lines,  Case  B. 

Note  that  the  advection  term  II,  blue  line,  is  always  negative  and  indicates  an  upward  flux  of  material. 
The  turbulent  diffusion  tenn.  III,  green  line,  can  be  of  either  sign,  with  negative  indicating  a  downward 
flux  of  material,  and  positive,  an  upward  flux.  The  growth  rate  tenn,  III,  black  line,  is  always  positive. 
The  three  tenns  must  balance  in  the  steady  state. 


In  Fig  3d,  the  no  turbulence,  pure  advection  (red  line)  regime,  where  RT  =  0  ,  advection  balances 

growth  rate.  For  the  turbulence  dominated  regime  of  Pe-  0  ,  Figure  3a,  turbulent  diffusion  balances 
growth  rate  for  both  case  A,  including  the  TIBI  term,  and  Case  B,  no  TIBI  tenn.  In  Figures  2  and  3a, 
where  Pe-  0  ,  note  the  constant  vertical  distribution  of  P  and  constant  vertical  distribution  of  terms 
III  and  IV  of  the  ADR  equation  (15).  This  results  from  vertically  uniform  mixing  in  the  steady  state. 
Turbulence  strongly  dominates  advection,  which  is  reflected  in  the  vertically  uniform  PDF,  F  =1. 
However,  this  result  only  depends  on  the  relative  value  of  the  distribution  of  terms  III  and  IV  of  the 
ADR  equation  (15).  This  results  from  vertically  uniform  mixing  in  the  steady  state.  Turbulence 
strongly  dominates  advection,  which  is  reflected  in  the  vertically  unifonn  PDF,  F  =1.  However,  this 
result  only  depends  on  the  relative  value  of  the  turbulent  time  scale  to  advective  time  scale  as  given  by 
the  definition  of  the  turbulent  Peclet  number, 


Pe  = 


i—r 

D 2 


a 


-i 


aD2 

K 


and  not  on  the  absolute  intensity  of  the  turbulent  field,  i.e.  value  of  ic  Ddone.  We  also  see  in  Figures 
3a  and  3b  for  turbulent  Peclet  numbers  Pe  -  0,1  that  for  Case  A,  TIBI  tenn  included,  (solid  green 
line),  at  all  depths,  turbulent  diffusion  fluxes  phytoplankton  toward  the  base  of  the  mixed  layer.  This 
downward  flux  of  phytoplankton  interacts  with  the  incoming  upward  moving  seed  nutrients  and  the 
growth  rate  (black  lines).  However  at  Pe  =  10  there  is  a  change  in  sign  of  turbulent  diffusion  with 
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depth  at  z  « .65  .  This  results  in  an  upward  turbulent  flux  of  phytoplankton  for  z  >  .65  .  The  upward 
flux  diminishes  the  growth  rate  in  that  depth  range,  as  indicated  in  Figure  3c.  Compare  in  Figure  3c  the 
black  solid  line  to  the  black  dashed  line,  the  fonner  being  with  TIBI,  the  latter  no  TIBI.  Thus,  the  TIBI 
tenn  affects  growth  rate  through  its  feed  back  with  the  turbulent  diffusion  tenn.  The  advection  term  is 
also  affected  when  the  TIBI  tenn  is  included.  Relative  to  not  including  TIBI,  it  remains  negative,  but 
upwardly  transporting  less  material  when  the  diffusion  tenn  III  is  positive,  and  more  material  when  the 
diffusion  tenn  is  negative.  These  effects  are  also  present  for  smaller  values  of  Pe,  where  turbulence 
dominates,  and  result  in  the  different  total  integrated  phytoplankton  growth  in  the  TIBI  versus  non 
TIBI  case.  See  Figure  2,  blue  and  green  solid  and  dashed  lines,  respectively. 

RESULTS 

In  the  late  1990s  Allan  Robinson  developed  a  theory  of  NPZ  interaction  in  a  laminar  upwelling  flow 
field.  His  approach  was  to  use  the  advection  reaction  equation  (AR)  and  obtain  solutions  for  the 
evolution  of  N,  P,  Z  in  a  Lagrangian  coordinate  system  (Robinson,  1997,  1999).  The  Robinson  theory 
was  extended  to  turbulent  flow  by  Goodman  and  Robinson  (2008)  by  using  a  probability  density 
function  (PDF)  on  the  solution  to  the  AR  equation.  The  PDF  employed  was  associated  with  a  random 
walk  undergoing  perfect  reflection  at  the  top  and  bottom  of  the  mixed  layer.  Bayes’  theorem  was  used 
to  express  the  PDF  in  Eulerian  coordinates. 

A  simple  bilinear  NP  turbulent  upwelling  interaction  case  was  examined  with  this  approach.  It  was 
shown  that  the  key  non  dimensional  parameter  describing  the  evolution  of  the  primary  production  was 
the  turbulent  Peclet  number,  the  ratio  of  the  advective  to  turbulent  time  scale.  It  was  also  observed  that 
the  Turbulence  Induced  Biodynamical  Interaction  (TIBI)  term  could  not,  in  general,  be  neglected.  The 
TIBI  term  arises  from  the  effect  of  turbulence  on  the  non  linear  part  of  biodynamical  interaction  and  is 
distinct  from  that  of  turbulent  mixing. 

However,  the  PDF  used  by  Goodman  and  Robinson  (2008)  did  not  satisfy  the  AR  equation  nor  did 
the  resulting  N,P,Z  satisfy  the  advective  diffusion  reaction  (ADR)  equation  .  This  resulted  in  an 
inability  of  that  approach  to  be  used  to  examine  the  role  of  the  TIBI  in  the  commonly  used  ADR 
models  of  NPZ  interaction. 

Allan  and  I,  just  prior  to  his  untimely  death,  developed  an  approach  ,  given  here  and  in  Goodman 
(2011)  to  obtain  a  PDF  which  does  satisfy  the  AD  equation  and  also  results  in  the  mean  biodynamical 
state  variables  N,P,Z  satisfying  the  ADR  equation.  This  approach  was  applied  to  the  bilinear  NP 
interaction  example  considered  in  previous  manuscripts  (Robinson  (1999)  and  Goodman  and  Robinson 
(2008),  with  particular  emphasis  on  the  role  of  the  TIBI  tenn.  Except  for  extremely  high  Peclet 
numbers,  Pe  »  10,  which  conespond  to  advection  effects  dominating  that  of  turbulence,  the  TIBI 

term  results  in  a  very  significant  contribution  to  the  mean  phytoplankton  profile  P  .  Neglect  of  the 

TIBI  tenn,  as  seen  in  figure  2,  results  in  an  overestimate  of  P  with  the  overestimate  increasing  with 
decreasing  Pe.  Not  including  the  TIBI  term  in  the  ADR  equation  also  greatly  alters  the  overall  role 
that  turbulence  plays  in  determining  the  mean  phytoplankton  profile  and  in  contributing  to  total 

phytoplankton  production.  As  indicated  in  figure  2,  for  the  limiting  case  of  Pe  =  oo  ,  not  including 

-1  _ 

TIBI  tenn,  results  in  the  mean  total  phytoplankton  production  in  the  steady  state,  PT  =  J  dzP  ,  being 
overestimated  by  approximately  40  %.  Also,  as  shown  in  Figures  (2)  and  (3)  ,  the  vertical  distribution 
of  P  and  the  three  terms  of  the  steady  state  ADR  equation  (14)  are  altered  in  the  vertical  over  the 
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intermediate  range  of  Peclet  numbers  1  <  Pe  <10.  In  particular  the  turbulent  diffusion  tenn  (green 
lines,  solid,  no  turbulence  case,  dashed  lines,  turbulent  cases)  show  a  sign  reversal  at  depth  for 
\<Pe<  10. 

The  original  theoretical  framework  developed  by  Allan  Robinson  lead  to  a  PDF  approach  to  modeling 
the  TIBI  term  and  most  recently  to  a  PDF  which  satisfies  the  AR  equation  and  with  N,  P,  Z  satisfying 
the  ADR  equation.  This  allows  a  quantification  of  the  role  of  TIBI  in  the  ADR  equation  approach  and 
a  prescription  on  how  to  proceed  with  more  complicated  and  realistic  NPZ  models  in  turbulent  flows. 


<P> 


Figure  2  Steady  state  ( t  — >  oo,  —  =  0  )  vertical  profiles  of  <P>  for  the  Peclet  number, Pe,  indicated  in 

dt 

the  figure,  using  a  =.l. .  Case  A,  the  solid  color  coded  lines  are  the  solutions  for  <P>  obtained  by 
using  (27a)  with  (26b)  as  the  PDF,  TIBI  term  included.  Case  B,  the  dashed  color  coded  lines,  the 
solution  to  the  ADR  equation  (28)  with  the  TIBI  term  neglected,  R  =  R0. 
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Figure  3  Steady  state  (t—>  oo, —  =  0 )  contribution  of  the  terms  of  equation  (28)-  advection,  term 
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II,  blue  lines;  turbulent  diffusion,  term  III,  green  lines;  and  the  growth  rate  (reaction)  ,  R,  term 
IV,  black  lines.  As  in  figure  1,  the  solid  lines  refer  to  Case  A,  the  dashed  lines,  Case  B.  Note  that: 
the  advection  term  II,  blue  line,  is  always  negative  and  indicates  an  upward  flux  of  material;  the 
turbulent  diffusion  term,  III,  green  line,  can  be  either  sign  with  negative  a  downward  flux  of 
material  and  positive  an  upward  flux;  the  reaction  (growth  rate)  term,  III,  black  line,  is  always 
positive.  The  three  terms  must  balance  in  the  steady  state.  Four  cases,  each  of  different  Peclet 

number,  Pe,  are  shown. 
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